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1.  INTRODUCTION 


Acoustic  backscatter  from  the  seafloor  has  generated  considerable 
interest  over  the  past  several  decades.  The  importance  of  seabed  acoustics  has 
recently  grown  due  to  current  emphasis  on  mine  warfare  in  the  littoral 
environment.  Since  there  is  at  present  no  general  agreement  on  the  physics  of 
the  scattering  process,  a  lot  of  models  exist,  based  on  several  different 
hypothetical  scattering  mechanisms.  They  include  models  based  on  scattering  at 
the  sediment  interface,  like  those  of  Patterson1  and  Clay  and  Medwin,2  as  well  as 
the  volume  scattering  models  of  Nolle  et  al.,3  Stockhausen,4  and  Ivakin  and 
Lysanov.5  Jackson  et  al.6"8  developed  a  two-component  model  that  treats 
interface  roughness  scattering  and  volume  scattering  separately.  Applying  this 
model  to  several  sites  suggests  that  both  components  are  important,  each 
dominating  the  total  backscatter  prediction  in  half  of  six  cases. 

Jackson’s  model  uses  six  parameters  to  characterize  the  seabed,  two  of 
which  describe  the  interface  roughness  and  three  of  which  describe  acoustic 
propagation  through  the  sediment,  modeled  as  an  acoustic  fluid.  The  remaining 
parameter  is  an  empirically  determined  quantity  which  specifies  the  level  of 
volume  scattering  within  the  sediment.  No  attempt  is  made  to  specify  the  actual 
volume  scattering  mechanisms  involved. 

In  1990,  Chotiros9  began  a  comprehensive  compilation  of  all  available 
shallow  grazing  angle  backscatter  data  published  in  the  literature.  The  results 
appeared  to  show  some  statistically  significant  trends  that  suggest  specific 
volume  scattering  mechanisms,  which  might  be  modeled  physically.  Two 
hypothetical  scattering  mechanisms,  involving  scattering  from  sediment  grains 
and  from  trapped  bubbles,  appear  to  fit  the  general  trends  in  the  data.  At  present 
the  data  set  is  expanding  and  the  trends  remain. 

In  1992  the  results  of  three  acoustic  penetration  experiments  at  sea10,11 
and  in  the  laboratory12  became  available.  They  suggested  that  sandy  sediments 
might  best  be  modeled  structurally  via  the  Biot  poroelastic  theory,  particularly  for 
high  frequencies  and  shallow  grazing  angles. 
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In  1993  development  of  a  new  backscatter  model13  was  begun  with  the 
intention  of  taking  advantage  of  the  most  recent  experimental  information.  The 
result  is  ARLUT’s  current  bottom  backscatter  model  called  BOGGART.  It  is 
intended  for  minehunting  applications  and  therefore  optimized  for  shallow  grazing 
angle,  high  frequency  behavior.  Like  Jackson’s  model,  BOGGART  includes 
volume  and  interface  scattering  components.  It  is  different  in  its  modeling  of  the 
sediment’s  structure  and  in  the  scattering  mechanisms. 

BOGGART  is  unique  among  current  scattering  models  in  two  important 
ways.  First,  the  sediment  is  modeled  structurally  via  the  Biot  theory.  Second,  the 
volume  scattering  is  modeled  physically  rather  than  empirically.  BOGGART 
consists  of  three  modular  components  that  predict  partial  scattering  strengths 
due  to  sediment  grains,  gas  bubbles,  and  interface  roughness.  Hence  the  name, 
which  is  an  acronym  for  “bottom  grain  gas  and  roughness  technique”. 

This  report  is  intended  to  serve  as  a  user's  manual.  Section  2  is  an 
overview  of  the  theory.  In  Section  3,  BOGGART'S  structure  is  described. 
Section  4  is  a  user’s  description  of  the  model's  operation. 
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2.  OVERVIEW  OF  THEORY 


This  section  discusses  the  underlying  concepts  upon  which  BOGGART  is 
built.  As  stated  in  the  introduction,  BOGGART  differs  from  other  current 
backscatter  models  because  it  incorporates  the  Biot  theory  to  describe  the 
seabed’s  structure,  and  because  the  volume  scattering  mechanisms  are  modeled 
physically  rather  than  empirically.  Section  2.1  discusses  the  Biot  theory’s 
significance  in  modeling  backscatter.  Section  2.2  explains  BOGGART’S  physical 
scattering  mechanisms. 

2.1  BIOT  THEORY 

The  Biot  theory  differs  from  other  sediment  theories  in  that  it  treats  the 
sediment  as  a  two-phase  medium,  consisting  of  a  solid  skeletal  structure  of 
sediment  grains,  through  which  a  pore  fluid  flows.  There  exist  two  compressional 
waves,  one  consisting  of  pore  fluid  and  grains  moving  in  phase  and  a  second 
consisting  of  the  pore  fluid  and  grains  moving  against  one  another  out  of  phase. 
The  first  of  these  waves  propagates  fastest  and  is  called  the  Biot  fast  wave;  the 
other  is  called  the  Biot  slow  wave.  Slow  waves  generally  attenuate  rapidly  in 
comparison  with  fast  waves  and  are  therefore  considered  negligible.  In  general, 
therefore,  sediments  are  modeled  with  simpler  one-phase  models,  such  as 
elastic  or  fluid  models. 

Recent  experiments  suggest  that  there  are  important  special  cases  in 
which  slow  waves  are  not  negligible,  making  the  Biot  theory  most  applicable.  In 
these  experiments,  accomplished  both  at  sea  and  in  the  laboratory,  a  slow 
acoustic  wave  was  observed  in  sandy  sediment  with  sound  speeds  near 
1200  m/s.  The  wave  was  significant  at  frequencies  greater  than  10  kHz  and 
grazing  angles  less  than  20°.  The  existence  of  this  wave  cannot  be  explained  via 
fluid  or  elastic  theories,  but  can  be  modeled  as  a  Biot  slow  wave  with  a 
reasonable  set  of  input  parameters. 

Measurements  of  normal  incidence  reflection  coefficients  also  suggest  that 
sandy  sediment  is  structured  as  a  Biot  medium.  Figure  2.1  shows  several  such 
measurements,  alongside  theoretical  values  based  on  Biot  and  elastic  theory. 
The  Biot  model  appears  to  be  most  consistent  with  the  measurements. 
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Figure  2.1 

Comparison  of  published  normal  incidence  reflection 
measurements  with  viscoelastic  and  Biot  theories. 
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The  Biot  theory  used  by  BOGGART  was  originally  formulated  by  Stern, 
Bedford,  and  Millwater.14  Its  incorporation  in  the  model  is  described  in  Ref.  13. 

2.2  SEDIMENT  SCATTERING  MECHANISMS 

Figure  2.2  is  a  plot  of  all  shallow  grazing  backscatter  measurements 
available  in  the  literature.  This  collection  was  compiled  by  Chotiros9  and  includes 
data  from  the  USA,  Canada,  the  UK,  and  Russia,  beginning  with  some 
measurements  recorded  by  Urick  in  the  1950s.  Each  data  point  represents  a 
backscattering  strength  measurement  over  sediment  at  a  grazing  angle  of  10°. 
The  data  span  a  broad  range  between  an  upper  limit  defined  by  Lambert’s  rule 
and  energy  conservation,  and  a  lower  limit  defined  by  Nolle’s  experiments3  over 
finely  graded  and  degassed  laboratory  sands.  An  outstanding  feature  of  Nolle’s 
data  is  that  the  backscattering  strength  increases  steadily  with  normalized  grain 
size  (grain  diameter/wavelength). 

If  the  data  points  are  grouped  according  to  the  site  they  were  collected 
from,  two  trends  appear,  shown  in  Fig.  2.3.  In  this  figure,  data  points  from  a 
common  site  are  connected  with  lines.  In  group  1,  the  backscattering  strength 
increases  with  normalized  grain  size,  just  as  was  the  case  for  Nolle’s  data  over 
degassed  sands.  This  dependence  of  the  scattering  strength  on  the  size  of  the 
sediment  grain  suggests  that  the  grains  are  involved  in  the  scattering  process.  In 
group  2,  the  backscattering  strength  appears  to  have  a  relative  maximum  when 
the  grain  diameter  is  about  10'2S  wavelengths.  This  is  close  to  the  resonance 
diameter  of  gas  bubbles  in  water,  suggesting  that  trapped  bubbles  in  resonance 
might  be  a  dominant  mechanism. 

BOGGART’S  hypothesis  is  that  the  shallow  grazing  angle  scattering  data 
of  Figs.  2.2  and  2.3  are  dominated  by  two  scattering  mechanisms:  grain 
scattering,  which  dominates  the  data  of  group  1,  and  bubble  scattering,  which 
dominates  group  2.  At  larger  grazing  angles,  it  is  possible  that  interface 
roughness  might  become  important.  BOGGART  models  these  three  scattering 
mechanisms  in  separate  modules,  described  below. 
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Figure  2.2 

Experimental  measurements  of  bottom  backscattering  strength 
as  a  function  of  grain  size  at  a  grazing  angle  of  10°, 
from  all  published  sources. 
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Figure  2.3 

Experimental  measurements  of  bottom  backscattering  strength 
as  a  function  of  grain  size  at  a  grazing  angle 
of  10°,  groups  1  and  2. 
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2.2.1  Sediment  Grains 


A  theoretical  model  for  acoustic  scattering  from  sediment  grains  is 
currently  under  development.15  It  will  be  based  on  the  Biot  theory14,16'18  for 
acoustic  propagation  through  fluid  saturated  porous  media,  with  a  modification 
that  includes  the  effects  of  scattering.  In  lieu  of  this  scattering  model,  which  is 
not  complete,  BOGGART  models  sediment  grain  scattering  empirically,  based  on 
experimental  data  measured  by  Nolle  et  al.3  The  expression  for  backscattering 
strength  is 


BSg  =  30  min 


-0.4623,  log 


-5.6  +  20  log 


r  sine  ] 

iiiwj 


(2.1) 


where 

d  =  effective  grain  diameter  in  meters, 
k  =  acoustic  wave  number  in  water,  and 
0  =  grazing  angle. 

2.2.2  Trapped  Gas  Bubbles 

The  trapped  bubble  component  of  the  backscattering  strength  is  based  on 
resonance  scattering  from  bubbles  that  exist  in  a  fluid  saturated  porous  medium. 
Acoustic  propagation  within  the  medium  is  modeled  by  the  Biot  theory.  Bubble 
resonance  behavior  inside  sediment  pores  is  modeled  via  a  simple  modification 
to  existing  theory  for  bubbles  surrounded  by  unbounded  fluids.2,19,20  In  this 
modification,  the  pore  fluid  is  assigned  an  effective  density  that  differs  from  its 
actual  density  and  accounts  for  the  fluid’s  partial  confinement.  The  details  are 
described  in  Ref.  13. 

Once  the  bubble  scattering  strength  is  known,  an  expression  is  developed 
for  the  backscattering  strength  per  unit  volume  of  sediment  containing  bubbles. 
An  effective  sediment  interface  scattering  strength  is  obtained  by  assuming  the 
backscatter  from  the  interface  is  a  sum  of  contributions  from  all  volume  elements 
below.  The  resulting  expression  for  backscattering  strength  by  bubbles  is 
expressed  in  the  following  form. 
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BSb  =  1 0  log  x 
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Plt  +  Pls 
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dz 
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where 

obVf  =  sediment  volume  scattering  cross  section  -  fast  wave, 

Cbvs  =  sediment  volume  scattering  cross  section  -  slow  wave, 

Pine  =  incident  acoustic  pressure, 

Pit  =  refracted  acoustic  pressure  -  fast  wave, 

Pis  =  refracted  acoustic  pressure  -  slow  wave, 

Po  =  fluid  density  above  sediment, 

Pit  =  effective  pore  fluid  density  -  fast  wave, 

pis  =  effective  pore  fluid  density  -  slow  wave, 

c0  =  sound  speed  in  water  column, 

cn  =  fast  wave  sound  speed, 

c1s  =  slow  wave  sound  speed,  and 

z  =  depth  below  sediment  interface. 

2.2.3  Interface  Roughness 


(2.2) 


The  backscatter  from  interface  roughness  is  modeled  by  solving  the 
Helmholz-Kirchhoff  integral.  In  this  regard  it  is  similar  to  the  approach  used  by 
Jackson  et  al.6  Both  models  use  the  Kirchhoff  approximation  which  limits  the 
frequency  range  of  applicability.  The  applicable  range  has  been  estimated  by 
Chotiros  for  the  measured  roughness  spectra  of  ocean  sediments  by  Briggs.22 
However,  the  method  adopted  here  is  different  from  that  of  Jackson  et  al.  in  two 
important  respects.  (1)  A  numerical  approach  is  used  to  compute  the  scattering 
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strength  as  a  function  of  grazing  angle,  in  order  that  measured  wave  number 
spectra  of  roughness  from  real  ocean  sediments,  expressed  in  the  form  of  a 
structure  function,  may  be  used.  Jackson  et  al.,6  on  the  other  hand,  attempted  to 
fit  a  linear  power  law  to  measured  wave  number  spectra.  As  a  result,  their  ability 
to  model  was  limited  because  the  integral  would  only  converge  within  a  limited 
range  of  power  law  slopes.  (2)  The  Biot  theory  is  used  to  determine  the  sediment 
plane  wave  reflection  coefficient  whereas,  in  Jackson  et  al.,  the  sediment  is 
approximated  as  a  fluid.  When  compared  to  measured  values  of  reflection  loss 
from  ocean  sediments,  the  Biot  theory  has  been  shown  to  give  more  realistic 
results.  The  plane  wave  backscattering  strength,  which  is  a  simplified  version  of 
a  result  developed  by  Chotiros21  for  roughness  scattering  under  more  general 
conditions,  is  given  by  the  following  expression. 


BSr=  10  log 


where 

k  =  acoustic  wave  number, 

0  =  grazing  angle, 

Dz  =  structure  function  of  roughness,  and 

R  =  plane  wave  reflection  coefficient,  a  function  of  the  grazing  angle. 

2.2.4  Total  Backscattering  Strength 

The  total  backscattering  strength  is  the  incoherent  sum  of  the  above  three 
components. 

BS=10log^100'1BS8+100-1BSb+100'1BSrj  .  (2.4) 
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3.  RUNNING  BOGGART 


In  order  to  run,  the  BOGGART  application  must  be  accompanied  by  five 
input  files.  When  the  program  is  run,  it  will  produce  three  output  files.  Input  and 
output  files  are  described  in  detail  in  Sections  3.1  and  3.2. 

3.1  INPUT  FILES 


Input  parameters  for  BOGGART  are  entered  via  input  files  which  can  be 
edited  by  the  user.  A  file  called  "BOTTOM. INPUT"  governs  the  program.  An 
example  is  given  in  Fig.  3.1(a).  This  file  selects  four  input  files  from  which  input 
parameters  are  gathered,  as  well  as  a  user  defined  prefix  which  will  be  attached 
to  all  output  files.  In  the  example  shown,  all  output  file  names  will  start  with  the 
prefix  “test”.  The  four  input  files  selected  by  the  user  are  (1)  the  sonar 
parameters  file,  (2)  the  sediment  parameters  file,  (3)  the  structure  function  file, 
and  (4)  the  sediment  gas  file. 


An  example  sonar  parameters  file  is  shown  in  Fig.  3.1(b).  The  sonar 
frequencies  and  grazing  angles  are  selected  here.  The  quantity  NTIMES  is  the 
number  of  frequencies  to  be  selected.  They  will  have  logarithmic  spacing 
between  a  minimum  "FS4M"  and  a  maximum  "FS4X",  inclusive.  If  NTIMES 
equals  1,  then  the  modeled  frequency  will  be  the  minimum  ‘FS4M’. 

An  example  sediment  parameters  file  is  given  in  Fig.  3.2.  This  file 
contains  the  Biot  parameters  that  describe  the  sediment.  They  are  broken  down 
into  three  groups:  (1)  solid  parameters  that  describe  the  grain  material,  (2)  fluid 
parameters  that  describe  the  pore  fluid,  and  (3)  frame  parameters  that  describe 
the  skeletal  structure  of  the  grain  matrix.  Four  general  sediment  parameters  files 
have  been  created  to  model  different  sediment  types  as  follows. 


Sediment  Type 


Name  of  Sediment  Parameters  File 


coarse  sand 
medium  sand 
fine  sand 
mud 


sediment.sandlo 

sediment.sand20 

sediment.sand30 

sediment.mud80 
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#File  selections  for  BotBS3: 

#Enter  appropriate  file  names  for  input  parameters  and  output  results 

sonar .parameters . file _  sonar .parameters 

site . sediment . file _  sediment . sand30 

structure . function . f ile_  struct . sand9 . 5mmRMS 

sediment . gas .  file _  gas.  layers 

output .  file  .prefix _  test 

#There  will  be  3  out  put  files  with  the  above  prefix 
#NNN. BSf ile  will  contain  backscattering  strength  vs  grazing  angle 
#NNN. const  will  contain  a  copy  of  all  the  input  parameters 
#NNN.vel  will  contain  the  wave  velocities  in  the  sediment 


(a)  Example  of  input  file  "BOTTOM.INPUT". 


sonar  parameters 


'  NTIMES-  .NUMBER.  OF.  FREQUENCIES _ '  1 

•  FS4M-  .  MIN .  FREQUENCY .  <kHz>  . _ '  30. 

'FS4X-.  MAX.  FREQUENCY.  <kHz>. _ '  30. 

'  NTHETA .  -  .  NUMBER .  OF  .  GRAZING .  ANGLES _ '  10  . 

'  THETMN-  . MIN .  GRAZING .  ANGLE .  <DEG>  . _ '  0  . 

'  THETMX-  .  MAX .  GRAZ  ING .  ANGLE .  <DEG>  . _ '  90. 


(b)  Example  of  sonar  parameters  file. 


Figure  3.1 
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sediment  parameters:  sand  MGS=  30 

Chotiros,  'Biot  model  of  sandy  sediments'. 

submitted  to  JASA; 

Boyle '  s  dissertation 

solid  parameters 

' RHOS- . SEDIMENT . GRAIN . DENSITY . <kg .mA-3>_ ' 

2650. 

' XKR- . . SEDIMENT . GRAIN . BULK . MOD . <Pa> 

7  .  e9 

' MEAN . GRAIN . PHI <MPHIs> 

3.0 

’ STAND . DEV . PHI<SIGPHIs> 

1.0 

fluid  parameters 

'  RHOF-  .  UPPER .  FLUID .  DENSITY .  <ka  . mA-3> 

1000 

' XKFO - . UPPER . FLUID . BULK . MOD . <Pa> 

2 . 25e9 

' RHOFBLO- . PORE . LIQUID . DENSITY . <kg . mA -3  >_ ' 

1000 

' XKFOBLO- . PORE . LIQUID . BULK . MODULUS . <Pa>_ ' 

2 . 25e9 

’ ETA- . . PORE . LIQUID .VISCOSITY . <kg .mA-l . S> ' 

1.00e-3 

' WATER . DEPTH<M> 

60. 

frame  parameters 

' B- . . . . SEDIMENT . POROSITY . <-> 

.50 

' XMUO - . SEDIMENT . FRAME . SHEAR . MOD . . <Pa> 

. 2  61e+8 

' DMUO - . SEDIMENT . FRAME . SHEAR . LOG . DECR . < - > ' 

.15 

' XKBO - . SEDIMENT . FRAME . BULK . MOD . <Pa> 

5 . 3e9 

' DKBO - . SEDIMENT . FRAME . BULK . LOG . DECR 

.15 

'  Vbp 

40.0 

Figure  3.2 

Example  of  sediment  parameters  file. 
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A  user  also  has  the  option  of  creating  a  new  sediment  parameters  file,  if  needed, 
and  selecting  it  inside  the  BOTTOM.INPUT  file. 

Figure  3.3  is  an  example  structure  function  file.  In  this  file,  sediment 
roughness  parameters  are  specified.  Two  generic  input  files  exist  at  present, 
which  apply  to  different  roughness  scales: 

rms  Roughness  Name  of  Structure  Function  File 

1.1  mm  struct.sandl.lmmRMS 

9.5  mm  struct.sand9.5mmRMS 

The  user  has  the  option  of  creating  a  new  structure  function  file  if  necessary,  and 
selecting  it  in  the  BOTTOM.INPUT  file. 

Figure  3.4  is  an  example  of  a  sediment  gas  file.  The  first  four  entries 
describe  the  physical  properties  of  the  gas.  At  the  bottom  of  the  file  is  the 
layered  gas  fraction  profile.  Each  row  consists  of  the  depth  of  the  top  of  the 
layer,  followed  by  the  layer’s  gas  fraction.  The  last  layer  is  assumed  to  be  an 
infinite  half  space.  There  is  no  limit  to  the  number  of  layers  which  can  be 
specified. 

3.2  OUTPUT  FILES 

BOGGART  produces  three  ASCII  output  files.  The  filenames  consist  of  a 
prefix,  specified  by  the  user  in  the  BOTTOM.INPUT  file,  and  the  suffix  “.const”, 
“.BS”,  or  “.vel”. 

The  “nnn. const”  file  is  a  readback  of  BOGGART’S  input  parameters.  An 
example  is  given  in  Fig.  3.5.  The  purpose  of  this  file  is  to  provide  a  complete 
record  of  the  input  parameter  values. 

The  “nnn.BS”  file  is  a  table  of  computed  backscattering  strengths  as  a 
function  of  grazing  angle  and  frequency.  Figure  3.6(a)  is  an  example.  The  first 
four  lines  are  a  list  of  the  input  files.  The  next  two  are  a  header,  followed  by  the 
data,  which  are  arranged  in  columns.  The  first  column  is  the  acoustic  frequency 
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Sand:Briggs  OCEAN  ENG  1989  14(4) ,Mission  bay  II  (isotropic) 
without  volume  scatter 

26  !  number  of  structure  function  samples  in  look  up 

table 
0  0 

0.000329222. 5365E-09 
0.000493835  .6758E-09 
0.000740741 .2616E-08 
0.001111112  .7638E-08 
0. 001666675.  88E-08 
0.0025  1 . 1877E-07 

0.00375  2 . 2481E-07 

0.005625  4 . 1234E-07 

0.0084375  7 . 5914E-07 
0.012656251. 3 663 E- 06 
0.018984382  .3 644E- 06 
0.028476563  .8765E-06 
0.042714845. 8 97 E- 06 
0.064072278. 2 72 4E-06 
0.0961084  1 . 1138E-05 
0.1441626  1 . 5345E-05 
0.2162439  2 . 1896E-05 
0.324365843  .1934E-05 
0.486548774  .6889E-05 
0.729823156. 8 19 6E-05 
1.094734739. 6382E-05 
1.642102090.00012911 
2.463153130.00015898 
3.6947297  0.000176 
5.542094550.00017976 


Figure  3.3 

Example  of  structure  function  file. 
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gas  parameters:  air,  medium  density 


' JHcpg-SPEC. HEAT. OF. GAS. .CONST. PRESSURE_  '  0.629 

•JHrhog-GAS.  DENSITY.  <kg.mA-3> _  '  1.22 

'XKG- .  .GAS.  BULK.  MOD.  .<Pa> _  '  2.48e5 

' DEPTH . OF . TOP . OF . LAYER . (M) , GAS . FRACTION : ’ 

0 . 0,  IE-3 
0.1, IE-5 


Figure  3.4 

Example  of  sediment  gas  parameters  file. 
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j  description (NO  SLASHES, SPACES  OR  COMMAS) data  field 

j  File  .Name .  of  .  out  .put .  results _ test 

;  sonar  parameter  input  file _  sonar .parameters 

!  sediment  site  input  file _  sediment . sand3o 

]  sediment  structure  input  file _  struct . sand9 , 5mmRMS 

;  gas  parameter  input  file _  gas. layers 


solid  parameters 

RHOS-. SEDIMENT. GRAIN. DENSITY. (kg.mA-3)_  2650.00 

XKR- . . SEDIMENT . GRAIN . BULK . MOD . (Pa) _  7.000000E+09 

MEAN  .GRAIN.  PHI  (rphis  ) _  3.00000 

STAND. DEV.  PHI  (dphis)__ _  1.00000 

fluid  parameters 

RHOF-. UPPER. FLUID. DENSITY. (kg.mA-3) _ _  1000.00 

XKF0-.  UPPER.  FLUID.  BULK.  MOD.  (Pa) _  2.250000E+09 

RHOFBLO- .PORE. LIQUID. DENSITY. (kg.mA-3)_  1000.00 

XKF0BLO-. PORE. LI QUID. BULK. MODULUS. (Pa)_  2 .250000E+09 

ETA- . .PORE. LIQUID. VISCOSITY. (kg.mA-l.s)  1.000000E-03 
WATER .  DEPTH.. _  60.0000 

gas  parameters 

JHcpg- SPEC. HEAT. OF. GAS. . CONST . PRESSURE_  1000.00 

JHckg -HEAT.  CONDUCTIVITY.  GAS _ .62  9000 

JHrhog- GAS  .DENSITY.  (kg.mA-3) _  1.22000 

XKG-  .  . GAS  .  BULK . MOD .  .  (Pa) _  248000  . 

Vbp _  40.0000 

frame  parameters 

B-  .  .  .  .SEDIMENT.  POROSITY.  (-) _  .500000 

XMU0- . SEDIMENT . FRAME . SHEAR . MOD . . ( Pa ) _  2 . 610000E+07 

DMU0- .SEDIMENT. FRAME. SHEAR. LOG. DECR. (-)  .150000 

XKB0-. SEDIMENT. FRAME. BULK. MOD. (Pa) _  5.300000E+09 

DKB0- .SEDIMENT. FRAME. BULK. LOG. DECR. . (-)  .150000 


dependent  constants 

B- . . . .SEDIMENT. POROSITY. (-) _ 

VMASS-. .VIRTUAL. MASS. CONSTANT. (-) 

PERM-. . . PERMEABILITY (mA 2) _ 

aa- . . . . PORE . SIZE . PARAMETER . (m) _ 

CC- ....  CONSTITUENT . PARAMETER . C _ 

CH- ....  CONSTITUENT . PARAMETER . H _ 

CM-. . . .CONSTITUENT. PARAMETER. M _ 


_  .500000 

_  1.50000 

_  1 . 134527E-10 

_  6 . 736548E-05 

(1.313699E+09, 1.5832 02E+08) 
(5 . 648118E+09,-l. 687774 E+ 08) 
(5 .387022E+09, -1 . 499875E+08 ) 


independent  variables 

NTIMES -NUMBER .  OF  .  FREQUENCIES _  1 

FS4M-MINIMUM.  FREQUENCY  (kHz) _  30.0000 

FS4X-MAXIMUM.  FREQUENCY  (kHz) _  3  0.0000 

NTHETA.- NUMBER.  OF.  GRAZING.  ANGLES _  10 

THETMN . -MINIMUM . GRAZING . ANGLE . DEG _  .000000 

THETMX. -MAXI MUM. GRAZING. ANGLE. DEG _  90.0000 


Gassy  layers:  top  boundary  and  gas  fraction 

number  of  gassy  layers _  2 

.000000  1 .  OOOOOOE-03 

.100000  1.000000E-05 

Structure  function  (m,mA2) 

number  of  structure  function  samples _  26 


.000000 
3 . 292200E- 
4 . 938300E- 
7 . 407400E- 
1.111110E-1 
1.666670E- 
2 . 500000E 
3 . 750000E 
5 . 625000E- 
8 . 437500E 
1 . 2  65625E- 
1 . 898438E- 
2 . 847656E- 
4 . 271484E' 
6.407227E- 
9 . 610840E- 
.144163  1 

.216244 
.324366 
.486549 
.729823 
1.09473 
1.64210 
2.46315 
3.69473 
5.54209 


000000 

04  2 . 536500E-09 

04  5 . 675800E-09 

04  1 . 261600E-08 

03  2 . 763800E-08 

03  5 . 880000E-08 

03  1 . 187700E- 07 

03  2 . 248100E-07 

03  4 . 123400E-07 

03  7 . 591400E-07 

02  1 . 3  66300E-06 

02  2 . 3  64400E-06 

02  3 . 876500E-06 

02  5 . 897000E-06 

02  8 . 272400E-06 

02  1 . 113800E-05 

. 534500E-05 
. 189  600E-05 
.  193400E-05 
.  688900E-05 
. 819600E-05 
.  638200E-05 
.  291100E-04 
. 589800E-04 
. 7  60000E-04 
. 797  600E-04 


Figure  3.5 

Example  output  file  "test.const":  list  of  input  parameter  values. 
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sonar  parameter  input  file _  sonar .parameters 

sediment  site  input  file _  sediment . sand30 

sediment  structure  input  file _  struct . sand9 . 5mmKMS 

gas  parameter  input  file _  gas  .  layers 

FREQ  | G  ANGL | FORWARD |  BACKSCATTERING  STRENGTH (dB) 


kHz 

1  deg 

LOSS 

dB 

deg . 

TOTAL 

GAS 

GRAINS 

ROUGH  X 

30, 

.000 

.0 

.00 

-180. 

.00 

-100 

.00 

- 

INF 

- 

INF 

-100. 

,00 

ROUGHNESS 

30, 

.000 

10, 

.0 

-8, 

.49 

-123. 

,32 

-19, 

.15 

-19 

.15 

-74 

.27 

-100. 

,00 

GAS 

30, 

.000 

20, 

.0 

-11, 

.26 

-59. 

,81 

-13. 

.31 

-13 

.31 

-68 

.38 

-91. 

09 

GAS 

30, 

.000 

30, 

.0 

-10. 

.73 

-22. 

,37 

-11, 

.08 

-11 

o 

00 

-65 

.09 

-76. 

09 

GAS 

30. 

.000 

40. 

.0 

-11. 

.79 

-1. 

00 

-10, 

.36 

-10 

.36 

-62 

.90 

-50. 

40 

GAS 

30. 

.000 

50. 

.0 

-12. 

.68 

84 

-10. 

.11 

-10 

.11 

-61 

.38 

-45. 

15 

GAS 

30. 

.000 

60. 

.0 

-12. 

.86 

1. 

00 

-9. 

.81 

-9 

.82 

-60 

.31 

-37. 

84 

GAS 

30. 

.000 

70. 

.0 

-12. 

.90 

1. 

04 

-9. 

.60 

-9 

.64 

-59 

.61 

-29. 

69 

GAS 

30. 

.000 

80. 

.0 

-12  . 

.90 

1. 

05 

-9. 

.22 

-9 

.50 

-59 

.20 

-21. 

37 

GAS 

30. 

.000 

90. 

.0 

-12  . 

.89 

1. 

05 

11. 

,65 

-9 

.45 

-59 

.07 

11. 

62 

ROUGHNESS 

(a)  Example  output  file  "test.BS":  backscattering  strengths. 


at  ( 1 . 57080 000000 ) degrees  grazing  and  30.0000kHz: 


f  ast .  wavespeed  (VELlo) _  1877.46 

slow,  wavespeed  (VEL2o) _  1288.32 

shear  .wavespeed  (VELso) _  131.239 

fast .  absorption  (ABSRPlo) _  34.3029 

slow,  absorption  (ABSRP2o) _  59 . 1957 

shear .  absorption  (ABSRPSo) _  398.683 


(b)  Example  output  file  "test.vel":  predicted  sediment  sound  speeds 
and  absorptions  in  m/s  and  dB/m. 


Figure  3.6 
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in  kHz.  The  second  is  the  grazing  angle  in  degrees.  The  third  and  fourth 
columns  are  the  forward  loss  and  phase  angle  of  the  reflection  coefficient.  The 
next  four  columns  are  backscattering  strengths  in  dB;  the  fifth  column  is  the  total 
predicted  backscattering  strength,  followed  by  the  individual  backscattering 
strength  components  for  trapped  gas,  sediment  grain,  and  interface  roughness 
scattering.  In  the  ninth  column,  the  dominant  scattering  mechanism  is  specified: 
either  “gas”,  “grains”,  or  “roughness". 

The  “nnn.vel”  file,  shown  in  Fig.  3.6(b),  is  a  list  of  computed  sound  speeds 
and  attenuations  within  the  sediment.  They  are  based  on  the  Biot  theory  and 
therefore  include  values  for  fast,  slow,  and  shear  waves.  The  purpose  of  this  file 
is  to  provide  acoustic  wave  speeds  and  attenuations  within  the  sediment  which 
may  be  useful  in  understanding  the  paths  by  which  the  sound  energy  engages 
the  scattering  mechanisms. 


19 


This  page  intentionally  left  blank. 


20 


4.  BOGGART  MODEL  STRUCTURE 


BOGGART  is  written  in  FORTRAN.  It  consists  of  the  driver  BotBS6  and 
22  subroutines  as  illustrated  in  Fig.  4.1.  The  subprograms  perform  the  following 
functions. 

BotBS6:  Serves  as  the  driver  for  BOGGART.  BotBS6  opens  input  files 

and  calls  subroutines  DATREAD3,  GASREAD,  and 
STRCTREAD  to  read  them.  It  then  calls  subroutines  BIOTK, 
SANDBS,  SURFBS,  GRAINBS,  and  ROUGBS  to  perform 
computations  and  write  the  results  to  the  appropriate  output 
files. 


DATREAD3:  Reads  in  all  input  parameters,  except  those  describing  the 

trapped  gas  or  the  interface  roughness. 


GASREAD:  Reads  in  parameters  pertaining  to  the  trapped  gas,  such  as  its 

density  and  bulk  modulus. 

STRCTREAD:  Reads  in  parameters  that  describe  the  interface  roughness. 

BIOTK:  Calculates  wave  numbers,  sound  speeds,  and  attenuations 

associated  with  Biot  fast  and  slow  waves  in  the  sediment  pore 
fluid.  It  calls  four  other  subroutines:  GASMOD,  INITCON2, 
TXRCOF,  and  RHOEF. 


GASMOD:  This  routine  exists  for  the  future  incorporation  of  an  algorithm  to 

compute  the  bulk  modulus  of  a  gassy  fluid.  At  present  no 
computations  occur  in  GASMOD,  although  two  variables  used 
elsewhere  in  BOGGART  are  defined.  CKM  is  the  complex 
compressibility  of  the  bubbly  pore  fluid,  presently  set  equal  to 
the  real  compressibility  of  the  fluid  without  bubbles.  RHOM  is 
the  density  of  the  bubbly  fluid,  presently  set  equal  to  the 
compressibility  of  the  fluid  without  bubbles. 
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Figure  4.1 

BOGGART  program  structure. 
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INITC0N2: 

MMKELO: 

MMKELD: 

CONST: 

WVNOS: 

TXRCOF: 

CTRIANG: 

RHOEF: 

SEDPRES: 

SANDBS: 


This  routine  contains  the  Biot  model’s  computations.  It  calls 
three  other  subroutines:  MMKELO,  MMKELD,  and  CONST. 

Calculates  the  Kelvin  functions  ber(x)  and  bei(x)  required  in  Biot 
model. 

Calculates  the  Kelvin  function  derivatives  ber’(x)  and  bei’(x) 
required  in  Biot  model. 

Calculates  some  of  the  constants  which  are  carried  through  the 
Biot  equations. 

Calculates  vertical  components  of  the  complex  wave  numbers  of 
the  fast,  slow,  and  shear  waves  within  the  sediment.  They  are 
assigned  the  names  “CL1”,  “CL2”,  and  “CLS”,  respectively. 

Calculates  Biot  transmission  and  reflection  coefficients  by 
Gaussian  elimination  of  the  linear  equations  (51  )-(54)  in  Ref.  14. 
Calls  subroutine  CTRIANG. 

Triangulates  a  complex  matrix  of  dimension  MxM. 

Calculates  fast  and  slow  wave  effective  densities  of  the  pore 
fluid  by  solution  of  Eqs.  (37)  and  (38)  in  Ref.  13.  It  calls  the 
subroutine  SEDPRES. 

Calculates  sediment  normal  tractions  and  pore  fluid  pressures  at 
the  fluid/sediment  boundary. 

Finds  the  volume  backscattering  cross  sections  abvs  and  obVf, 
given  the  appropriate  pore  fluid  effective  density  calculated  by 
the  subroutine  RHOEF. 
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GRANDIST: 

PORCONV: 

ROTATE: 

GASPRM: 

BAKSCAT: 

SURFBS: 

SURFBS: 


Calculates  shape  of  bubble  size  distribution  from  measured 
grain  size  distribution.  The  result  is  entered  into  an  array  called 
“BUBPOP”.  This  array  is  not  normalized.  The  total  volume  of 
gas  subtended  by  BUBPOP,  per  volume  of  pore  fluid,  is 
assigned  the  name  INTGR.  Calls  subroutines  PORCONV, 
ROTATE. 

Calculates  pore  size  density  function  by  convolving  grain  size 
density  function  with  a  pore  size  density  function  of  a  dense 
random  packing  of  hard  spheres,  according  to  Eq.  (14)  in 
Ref.  13. 

Rotates  a  real  one-dimensional  array  of  “NARRAY”  data  points 
by  “NROT’ points. 

Normalizes  BUBPOP  to  obtain  a  bubble  size  density  function 
called  BPNORM,  whose  total  gas  volume  equals  the  sediment 
porosity. 

Finds  the  backscattering  strength  per  unit  volume  of  sediment. 
BAKSCAT  is  called  twice  inside  SAN  DBS,  once  for  the  Biot  fast 
wave  and  once  for  the  slow  wave. 

Finds  effective  backscattering  strength  of  sediment  interface  by 
summing  contributions  from  the  sediment  volume  below.  The 
sediment  volume  is  broken  up  into  depth  cells,  from  which 
contributions  are  summed. 

Adjusts  the  depth  increment  of  the  cells  so  that  the  Biot  wave 
with  the  fastest  decay  rate  decreases  its  pressure  amplitude  by 
10%  from  the  top  to  the  bottom  of  a  single  cell.  The  integration 
is  carried  to  a  depth  determined  by  the  pressure  amplitude  of 
the  wave  with  the  slowest  decay  rate.  When  this  wave  has  1 0% 
of  its  amplitude  at  the  interface,  the  integration  is  stopped.  If  a 
medium  with  multiple  layers  is  modeled,  the  integration  is  not 
allowed  to  stop  until  the  last  layer  is  reached. 
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SURFBS: 

GRAINBS: 

ROUGBS: 


Calls  one  other  subroutine,  “SEDPRES”  to  get  the  refracted 
acoustic  pressures  at  the  water/sediment  interface. 

Calculates  the  backscattering  strength  contribution  due  to 
scattering  from  sediment  grains.  Empirical  modeling,  based  on 
the  experiments  of  Nolle  et  al.,3  is  currently  used  in  lieu  of  a 
theoretical  model,  which  is  under  development. 

Calculates  the  backscatter  component  due  to  interface 
roughness  by  evaluation  of  the  Helmholz-Kirchhof  integral. 
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